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Abstract 

Offline computation is an essential component in most multiscale model reduction techniques. How¬ 
ever, there are multiscale problems in which offline procedure is insufficient to give accurate representa¬ 
tions of solutions, due to the fact that offline computations are typically performed locally and global 
information is missing in these offline information. To tackle this difficulty, we develop an online local 
adaptivity technique for local multiscale model reduction problems. We design new online basis func¬ 
tions within Discontinuous Galerkin method based on local residuals and some optimally estimates. The 
resulting basis functions are able to capture the solution efficiently and accurately, and are added to 
the approximation iteratively. Moreover, we show that the iterative procedure is convergent with a rate 
independent of physical scales if the initial space is chosen carefully. Our analysis also gives a guideline 
on how to choose the initial space. We present some numerical examples to show the performance of the 
proposed method. 


1 Introduction 

In this paper, we develop an online local adaptivity technique for a class of multiscale model reduction 
problems. Many realistic applications involve solving problems that contain multiple scales and high contrast. 
Direct solution methods for these problems require fine-grid discretizations and result in large discrete systems 
that are computationally intractable. Common model reduction techniques perform the discretization of the 
problems on a coarse grid, which is much larger than the scales under consideration, with the aim of getting 
more efficient solution strategies. There are a variety of multiscale model reduction techniques based on 
numerical upscaling (e.g, mm) or multiscale methods (e.g., [H [5j [15j HU [22l [23l [25l [6j [TlJ 0 [9]). Most 
of the existing techniques are based on the so called offline construction. In particular, reduced models are 
computed in a pre-processing step, called offline stage, before the actual simulations, called online stage, are 
performed. For instances, some effective media are pre-computed for methods based on numerical upscaling 
and some multiscale basis functions are pre-computed for multiscale finite element methods. While these 
methods are effective in a wide variety of applications, there are still situations for which these methods 
are inadequate to give reliable solutions unless a large dimensional offline space is employed. Some of these 
situations involve external source effects and distant effects, which are ignored by most multiscale model 
reduction methods since they are typically based on local constructions. Therefore, it is evident that offline 
procedures are sometimes not enough to give efficient reduced models. Hence, it is the purpose of this paper to 
design a novel multiscale model reduction method. Our proposed method is based on a combination of offline 
technique and an online enrichment technique. The online technique is able to produce a reduced model 
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taking care of external sources and distant effects, without using global models. The online construction 
is also performed locally and adaptively in regions with more heterogeneities, giving very efficient reduced 
models. 

Our proposed method follows the overall idea of the Generalized Multiscale Finite Element Method 
(GMsFEM), which is introduced in T7] and is a generalization of the classical multiscale finite element 
method ([2Sj) in the way that the coarse spaces are systematically enriched, taking into account small scale 
information and complex input spaces. Instead of conforming finite element spaces as in QUO!!, we will 
use in this paper discontinuous Galerkin finite element spaces, which have some essential advantages (see 
use]) in multiscale simulations because it allows coupling discontinuous basis functions. The discretization 
starts with a coarse grid and a space of snapshot functions, which are defined on coarse elements. A space 
reduction is then performed to obtain a much smaller offline space by means of spectral decomposition. 
The spectral decomposition is performed locally on coarse elements, thus the functions in the offline space 
are in general discontinuous across coarse edges. The offline space is used as the approximation space for 
the interior penalty discontinuous Galerkin (IPDG) discretization on the coarse grid for the problem under 
consideration, giving our generalized multiscale discontinuous Galerkin method (GMsDGM). We remark 
that the offline space is computed only once in the pre-processing offline stage, and the same set of basis 
functions is used for any given source terms and boundary conditions. A-priori error estimate can be derived 
as in mm® showing that the error is inverse proportional to the first eigenvalue corresponding to the 
first eigenfunction that is not used in the construction of the reduced space. Since the aim of the paper is 
the new online locally adaptive procedure and its convergence, we will not discuss a-priori error estimate in 
this paper. 

The previous paragraph discusses the offline component of our method. As we discussed before, some 
new basis functions are necessary to capture certain behavior of the solution which cannot be captured by 
offline basis functions. For example, the solution may contain heterogeneities due to some distant effects 
and source terms, and these cannot be incorporated efficiently by offline basis functions before the solution 
is computed. Hence, it is the purpose of this paper to develop a technique to find new basis functions in 
the online stage. Our method consists of an iterative procedure. Given an approximate solution, some local 
residuals on coarse elements can be computed to reflect the amount of error in these coarse elements. These 
local residuals serve as indicators to locate regions, where new online basis functions are necessary. We will 
show that the projection of these residuals to the fine-grid can be used as new basis functions and that 
the energy-norm error has the most decay in a certain sense when these residual-based basis functions are 
included in the next solution process. In addition, we will show that this iterative procedure is convergent 
with a convergence rate independent of scales and contrast. In our analysis of convergence, we will show 
that it is essential to choose the appropriate space to begin the iterative procedure (cf. [16] [24] [19]). This 
initial space is computed in the offline stage and is obtained from a carefully design spectral problem. With 
this choice of the initial space, we show that one can obtain a very fast decay of errors by adding our online 
basis functions. Finally, we remark that there are offline adaptive enrichment strategies in the context of 
GMsFEM. In particular, in im offline adaptive procedure is developed and its convergence is analyzed 
using techniques in mug. This is an efficient method to adaptively enrich the offline space and is desirable 
for problems, where offline basis functions are good enough to capture the solution. On the other hand, we 
remark that other adaptive methods are available Also, we remark that online basis 

functions within continuous Galerkin GMsFEM is studied in [5j- 

The rest of the paper is organized in the following way. In the next section, we present the basic idea of 
GMsDGM and our online locally adaptive procedure. The method is then detailed and analyzed in Section 
[3] In Section [4] numerical results are illustrated to test the performance of this adaptive algorithm. Finally 
a conclusion is given in Section [5] 
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2 Preliminaries 

In this paper, we consider the following high-contrast flow problem 

— div(/t(a;) Vit) = / in D, (1) 

subject to the homogeneous Dirichlet boundary condition u = g on dD, where D is the computational 
domain and f(x) is a given source term. We assume that the coefficient k(x) is highly heterogeneous with 
very high contrast. For the convenience of our analysis, we also assume that k(x) is bounded below, that 
is, k(x) > 1. Due to the heterogeneity and high contrast of n(x), very fine meshes are necessary to obtain 
accurate numerical solutions. It is therefore crucial to develop a numerical scheme with a low dimensional 
approximation space for the efficient approximation of (HD- 

Next, we present some notations needed for the constructions of our scheme. Consider a given triangu¬ 
lation T h of the domain D with mesh size H > 0. For convenience, we assume the domain D is rectangular 
and that the triangulation T H consists of rectangles. We call T H the coarse grid and H the coarse mesh 
size. Elements of T H are called coarse grid blocks and we use N to denote the number of coarse grid blocks. 
The set of all coarse grid edges is denoted by £ H . See Figure [Q for an illustration. We also introduce a finer 
triangulation T h of the computational domain D, obtained by a conforming refinement of the coarse grid 
T h . We call T h the fine grid and h > 0 the fine mesh size. 
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Figure 1: Illustration of a coarse neighborhood and a coarse element. 


Now we present the framework of our GMsDGM. The methodology consists of two main ingredients, 
namely, the construction of local basis functions and the global coarse grid level coupling. For the coarse 
grid level coupling, we will apply the interior penalty discontinuous Galerkin (IPDG) method [3D]. We 
remark that other discretizations can also be used. Assume that Vr is a given approximation space defined 
on the coarse grid T H . Functions in Vr are piecewise polynomials within coarse grid blocks and are in 
general discontinuous across coarse grid edges. Following standard procedures, the IPDG method reads: 
find ur G Vr such that 

odg(ur,v) = (f,v), VvGVr, (2) 


where the bilinear form odg is defined as 


odg(«,«) = a H (u,v) - V) / ({kVu • n B }[u] + {kVu • + V) 7 / 

Ee£» jE E^ hjE 


(3) 
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(4) 


a H {u,v) 


Y a%(u,v), 
KgTh 



k,X7u • Vw, 


where 7 > 0 is a penalty parameter, he is a fixed unit normal vector defined on the coarse edge E £ S H . 
Note that, in (|3jl. the average and the jump operators are defined in the classical way. Specifically, consider 
an interior coarse edge E £ £ H and let K + and K~ be the two coarse grid blocks sharing the edge E. For 
a piecewise smooth function G with respect to the coarse grid T H , we define 


{G} = l(G+ + G~), [G] = G + — G~, on E, 


where G + = G \x+ and G~ = G\x- and we assume that the normal vector he is pointing from K + to K~. 
Moreover, on the edge E, we define k = (kk+ + k k-)/ 2, where Kk± is the maximum value of k over AT . 
For a coarse edge E lying on the boundary dD, we define 


{G} = [G] = G, and n = kk on E, 


where we always assume that tie is pointing outside of D. 
For our analysis, we define the DG-norm as 


IdG ““ a H(u, u) + 


V 1 


c[u] 2 . 


EeE 1 


Then, the following continuity and coercivity of the bilinear form odg hold. For completeness, we include a 
proof of this result in the Appendix. 

Lemma 2.1. Assume that the penalty parameter 7 is chosen so that 7 > G 2 1V . The bilinear form odg 
defined in m is continuous and coercive with respect to the DG-norm, that is, 


a DG (u,v ) < ai ||u||dg IMIdg, (5) 

a DG {u,u) > aolMlijc ( 6 ) 

for all u,v £ Vh, where ao = 1 — Gi nv 7 - ^ > 0 and ai = 1 + G; n v 7 _ ^. 

One main result of the paper is a convergence estimate of an adaptive procedure for the problem ([2]). For 
this purpose, we will compare the multiscale solution uh to a fine-scale solution Uh defined in the following 
way. We first let 

v£a = {v£ L 2 (D) : v\ K £ V h {K)}, 

where V h (K) is the space of continuous piecewise bilinear functions defined on K with respect to the fine 
grid. The fine-scale solution Uh £ V^g defined as the solution of the following 

odgK.w) = (f,v), Vv£V£ g . (7) 


It is well-known that Uh gives a good approximation to the exact solution u up to a coarse grid discretization 
error. 

The second main component of our method is the construction of local basis functions, which contains 
two stages, namely the offline stage and the online stage. In the offline stage, a snapshot space I/ l > sna P is 
first constructed for each coarse grid block Ki £ T H . The snapshot space contains a rich space of basis 
functions, which can be used to approximate the fine-scale solution defined 0 with a good accuracy. A 
spectral problem is then solved in the snapshot space y** sna P and eigenfunctions corresponding to dominant 
modes are used as the basis functions. The resulting space is called the local offline space V l ’ oS for the i-th 
coarse grid block Ki. The global offline space V oS is then defined as the linear span of all these V l ' oS , for 
i = 1,2, • • • ,N. This global offline space will be used as the initial space of our method. We denote 
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this initial space as Vf 1 / 1 . Using the initial space, an initial solution u$ can be computed by solving (|2jl. 
Local residuals in coarse grid blocks can then be computed based on the initial solution u^P. In coarse grid 
blocks with large residuals, new basis functions are computed and added to the approximation space. This 
procedure is continued until certain tolerance is reached. Next, we present a general outline of the method. 

Assume that the initial space is given and the initial solution is computed. For any m > 0, we 
repeat the following until the solution satisfies certain tolerance requirement. 

Step 1 : Solve © using the space V^ n> to obtain the solution u^ £ V^ n \ 

Step 2: Compute local residuals based on the solution 

Step 3: Construct new basis functions in regions, where the residuals are large. 

Step 4: Add these basis functions to V^"' > to form a new space y( m+1 ). 

In the following, we will give the details of Step 2 and Step 3. We will also explain how one chooses the 
initial space vj^\ 

3 Locally online adaptivity 

In this section, we will give details of our locally online adaptivity for the problem ©. As presented in the 
general outline of the method from the previous section, our adaptivity idea contains the choice of initial 
space as well as construction of new local multiscale basis functions. In the following, we will give the 
construction of these in detail. 

3.1 Initial space 

We present the definition of the initial space Let Xi be a node in the coarse grid T H , referred to as 

the i-th coarse node, for i = 1,2, ■ • • , N c> where N c is the number of nodes in the coarse grid T H ■ We will 
then define the i-th coarse neighbourhood as the union of all coarse grid blocks having the node Xi, see 
Figure Q] Moreover, for each coarse grid block K £ T H , we let xfj)i 3 = 1,2, 3,4, be the partition of unity 
functions, having value 1 at one vertex yj and value 0 at the remaining three vertices, where yj, j = 1 , 2 ,3,4, 
are the four vertices of K. Note that there is exactly one value of j such that the vertex y 7 is the same as the 
vertex Xi . In the case, we write x?\ = xf ■ One can use the standard multiscale basis functions or bilinear 
functions as the partition of unity functions. Note that we do not require any continuity of these partition 
of unity functions across coarse grid edges. The partition of unity functions are all supported on coarse grid 
blocks. Furthermore, we define the space V h (u>i ) by 

V\ui) = {v£ L 2 {iOi) : v\ K £ V h (K), K £ T H , K C W< }. 

That is, functions in V h (ui) are supported in Ui and belong to the space V h (K) for each coarse grid block 
I\ C uii- Note that there is no continuity condition across boundaries of coarse grid blocks. We consider 
V h (u>i) as the snapshot space in u>i, that is I/ z - sna P = V h (u ;*), and perform a dimension reduction through a 
spectral problem. For this purpose, we define £[ r be the set of coarse grid edges lying in the interior of 0 J t , 
and the following bilinear form 

a Ui (u,v)= a^(u,v)+ ^2 t f «HW, Vu,v £ V h (wi). (8) 

KgT h ,KCuh EGSf 1 JE 

Based on our analysis to be presented next, we solve the following spectral problem 

auH{u,v) = Xs Ui (u,v), Vv£V h (u l ), (9) 
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where 


( 10 ) 


s Ui (u,v)= / ^xf\ 2 uv + J «[xf] 2 MW, Vu,v eV h (u)i). 

kgt h ,kcuh^ k Eesf 1 Je 

We use the notations and to denote the fc-th eigenvalue and the fc-th eigenvector of the above 
spectral problem ©. Each eigenfunction corresponds to a function in ?/>£* G V h (ux) defined by 

Ui 

C = DW. 


where rit is the dimension of V h (ui ) and {wj*}™ 1 1 is a basis for V h (uji). In the above definition, is 

the j-th component of the eigenvector dAb 

For each Wj, we solve the spectral problem © and the first eigenfunctions are used to form the initial 
space. Each eigenfunction ip 1 ?* will be first multiplied by the partition of unity function xfx for each K C Wj, 
and is then decoupled across coarse grid edges to form 4 basis functions. In particular, the 4 new basis 
functions have support in one of the coarse grid block forming u>i and are zero in the other three coarse 
grid blocks forming Wj. For example, if K C u 4 , the basis function is xf'V’fe*- We W 1 'it e W ,off as the space 
spanned by > f° r a U ^ C w*. See Figure [2] for an illustration. The initial space vjf 1 is obtained by 

the linear span of all functions constructed in the above procedure. 



Figure 2: Illustration of the initial basis construction. Left: An eigenfunction ?/;£’" is defined in w*. Right: 
4 basis functions are obtained by splitting xf^fc 4 into 4 pieces, and each has support in K C Wj. 


3.2 Construction of online basis functions 

In this section, we will discuss the construction of our local online basis functions. The purpose is to add basis 
functions locally in some coarse neighborhoods to obtain rapidly decaying errors. Assume that the space 
at the m-th iteration and the corresponding solution u^ 1 ' 1 are given. For each coarse neighborhood w*, 
we define the local residual by 


Ri m \ v ) = l/» - odg( w^ n) ,-y), v G 


( 11 ) 
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where V 0 h (uJi) C V h (iOi) contains functions that are zero on du>i. The residual i?K can be seen as a linear 
functional defined on V 0 h (cui) with norm ||i?K|| defined by 


I*r>i- s»p 

vGV 0 h (uii) ll w ll^i 


and ||t;|| 2 . = a Ui {v,v). We will then find the new online basis function <f) G VgKi) by solving 


a Ui ((j),v) = Vu G V 0 h (u;i). 


( 12 ) 


The new online basis function ^ is added to V^" 1 to form V^ n+1 ' > . 

The motivation of finding a new basis function <j) by solving (fl2l) can be explained as follows. We define 
the A-norm by 

IMIa = o-dg(u, u), Vu G Vh. 

We notice that the A-norm is equivalent to the DG-norm ||it||DG by Lemma 12.II From ([2]) and J7]), we have 
the following Galerkin orthogonality condition 


odgK ~ v) = 0, Vu G Vjj m \ 

Thus, we see that the following optimal error bound holds 

W'U'h II A. — W^h ^ 11 A i VU £ Vj[ 

Notice that m and m hold for any m > 0. 

We will enrich the space V^"' 1 by adding a basis function <j) in the space V^(u!i) to form Vj , j " +l> 
(fHl) implies 


K -U { h +1) \\ 2 a < \\u h -u\\ 2 A , 


VuG f4 m+1) - 


Taking u = u\ r f nl + a<j>, for some scalar a , in (fl5l) . we have 

IK - u ( ™ +1) ||i < || u h - - Klli 


(13) 


(14) 


. First, 
(15) 


which implies 

II u h - u ( h +1) ||i < || u h - u^ !) || i - 2aa DG (u h - ,$) + a 2 \\(/)\\ 2 A . 

Taking a = odg(w/ 1 — </>)/||^>||ii we obtain 


/ (m) ,\2 

II u h - u^ +l) Hi < IK - vg> Hi - QdgK ||~|| M / ■ 


By the definition of the residual in (fill) , we see that (1161) becomes 


(16) 


^-4 m+1) lli<IK-K m ili- 


(m ) l|2 (C ^)) 2 


IKi 


(17) 


From m , we see that the quantity (i?K(^)) 2 /||0||i measures the amount of reduction in error when the 
basis function <j) is added in vjf 1 ^ to form vj^" +1 \ We will construct the function <j> G Vq K ) to obtain the 
most reduction in error. Thus, we find <j> G V 0 h (uii ) that maximizes R[ m \<p)/\\(j)\\A- Equivalently, we find 
cj> G V 0 h (uJi ) by solving 

v) = R\ m \v), Vu G V 0 h (cji). 

Notice that, we have used the fact that ||</>||a = || <j)\\ bJi when cj> G V 0 h (cOi). 
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3.3 Convergence of the adaptive procedure 

In this section, we analyze the convergence of the above online enrichment procedure. We begin our analysis 
at the inequality O- Notice that, this inequality can be written as 

ik - »sr +1, iii < ik - <4"’e - iiR.! m, ii 2 , IK 

when the basis function <f> is obtained as in m- 

On the other hand, we will show that the error ||— u^ 1) \\a can be controlled by the residual norm 
To do so, we consider an arbitrary function v £ V GG . Let i\ £ V h {uji) be the restriction of v in 04 , 
and let £ U* ,off be the component of w, in the offline space U I,off . By the GMsDGM ([2]), the fine-grid 
problem 0 and the Galerkin orthogonality (1131) . we have 

odg(m/i - ^\v) = a DG (u h - u^\v- u (0) ), Vu (0) G , 

where we define e and use the fact that U^ 0) C vjf 1 '* for all m > 0. By 0, we have 

odg(m/i - u'h^v) = (f,v - u (0) ) - a DG (u^ l \v - i/ 0) ). 

Using the property xfj) = 1 f° r a ll K G T H , 

av G (uh -u { h\v) = J2 {(f’X?j)(v - v^)) - a BG (u { x\x?j)(v ~ 

KeT H i=i 


Writing the above sum over coarse neighborhoods w*, we have 

N c 

«dgK - u{ h\v) = 

i= 1 KduJi 

For each coarse neighborhood cjj, we define the following modified local residual by 

Ri m \ v ) = ((/’ -« DG ( u iT ) iXf w ))> v£V h (uj i ). (19) 

KCuii 


The modified residual R^ n ' > 
in the following way 


can be seen as a linear functional defined on V h (uii) with norm 

||S(" 0 ||_ „ lm \Rt\v)\ 

II II SUp if 11 

v£V h (uji) || Xi v \\ui 


defined 


In the above definitions, Xi is considered to be defined only on K , and has zero value outside K. 
Using the definition of the modified residual R^ m \ we have 


N c 


oidgK - ll-Ri m) || II X] xf(v-^ 0) )\\ A 


( 20 ) 


i— 1 


KCoji 


where we used the fact that ]> Zkcoj X? ( v — v i°' > ) is zer0 on cC;. Using Lemma |2~T 


X?(v - W )IU < «i 2 II X?{v - W )I|dg- 




KtZuii 


( 21 ) 



By the definition of the DG-norm, 


K <•• k 0) )1 2 . 


llxf(«-«F)llDG = E [ + [nxf(v- 

KCUi JK e , ' e 

For each K C Wj, we have 

/ «|V(xf (w-v- 0) ))| 2 < 2 f «x?|V(i;-u,< 0) )| 2 + 2 [ «|Vxf|>-<4 0) ) 2 . 

JK JK JK 

For each e E £/*, we have 

/ «[*?(« - ^ 0) )F < 2 f nxfVlv vi 0) r + 2 f Ax?f{v - ^ (0) } 2 - 

J e J e J e 

Combining inequalities (l23l) and (l2dl) in (l22l) . we have 

llxf(^-'y|° ) )llDG<2||'y-wf ) ||i. +2( f 75|Vx< r | 2 (w — 0) ) 2 + tY 1 j^bcffb-vf} 2 ) (25) 

Jwi ee£f je 

where ||f||^. = a Wi (u,u). Using the spectral problem (0, we have 

Ik - f | 0) ||Ai = - Wi 0 ) ,u* - f| 0) ) < a Ui (vi,Vi) = Ikll^. 


and 


[ «|v%f l 2 k -^ 0) ) 2 + t 5Z [ K l^xff{v-v ( f ) } 2 = s Ui {v i -v i f\v i -v { ° ) ) < tid;— IMU- 

J Ui n ed£ aJe A Li + l 


( 22 ) 


(23) 


(24) 


Thus, (1251) and (1211) impleis 


Hence, EOl) becomes 


Wxf {v ~ vf ] )\\ 2 A <2ai(l + )|| 


Li+1' 


N c 

sdg(«/. - < ( y^2ai(l 

i=l 


"'Li + 1 


I R 


(m),|2\ 2 


N c i 

(EiHii) 1 . 


We remark that the above inequality holds for any v £ V^ G . Taking v = Vh — v^ 1 ' and using Lemma \2. II 
we finally obtain 

N e 

IK - u™ \\ 2 A < 2a 0 - 1 a 1 C 0 ]T (l + ||^ m) || 2 , (26) 

i= 1 Li + 1 

where Co = u\ax.x^T H n K and Uk is the number of vertices of the coarse grid block K. 

We define 

e=\\R { f Ci f/n 2 , and T ] 2 = 2a ( C 1 aiC 0 ^ (l + )||^-”' , || 2 . (27) 

i— 1 Li +1 

From (fl8l) and (l26l) . we see that the following convergence holds 

\\u h -u { H +1) \\ 2 A < (1 - 0 )\\u h - u^Wa- 


We summarize the above results in the following theorem. 
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Theorem 3.1. Let Uh be the solution of and u^"' 1 , m > 0, be the solution of (0|) in the m-th iteration. 
Then the following residual bound holds 


N c 

IK -«?>£ < 2%'a.C„E (l + )ll^ rn> ll 2 * 


(28) 


ii+1' 


Moreover, the following convergence holds 


\\ u h - u { h +1] \\a < (1 -0)\\u h -u { H ] \\ 


H IIA 


(29) 


where 6 is defined in \2 7| ). 

We remark that one can derive a priori error estimate for the error — u^”*||dGi for every m > 0. 
Since the purpose of this paper is an a posteriori error estimate (I2SD and the convergence of an adaptive 
enrichment algorithm (1291) . we will not derive a priori error estimate. 

Finally, we remark that by using more basis functions in the initial space V ^, the values of the eigenvalues 
\ff +l are larger. Thus, the value of 6 is further away from zero, and this fact enhances the convergence rate. 
In particular, the convergence rate is affected by the quantity A m ; n = mini<j<jv c The convergence is 

slow when A m i n is small (cf. 0 [24]). We also remark that one can add online basis functions in multiple 
coarse neighborhoods to speed up the convergence. Let S be the index set for which online basis functions 
are added in w* for i £ S. By using similar arguments as above, we obtain 

\\u h -u { ™ +1) \\ 2 A < (i-e)\\u h -u { ™ ] \\\ 


where 

0 = EK (m) IW- 

ies 


4 Numerical Results 

In this section, wc will present some numerical examples to show the performance of the proposed method. 
The implementation procedure of online adaptive GMsDGM is described below. First, we choose a fixed 
number of functions for every coarse neighborhood by solving the local spectral problem. This fixed number 
for every coarse neighborhood is called the number of initial basis. After that, we split these functions into 
the basis functions of the offline space such that each basis function is supported in one coarse grid block. 
We denote this offline space as V oS and set = V oS . 

The coarse neighborhoods are denoted by oJij, where i = 1, 2, • • ■ , N x and j = 1,2, • ■ ■ , N y and N x and 
N y are the number of coarse nodes in the x and y directions respectively. We consider I x ,odd and I x ,eve n 
as the set of odd and even indices from {1,2,- ,N X }. Similarly, I yt0 dd and Iy^ven are the set of odd 
and even indices from (1, 2, ■ • • , N y }. In each iteration of our online adaptive GMsDGM, we will perform 
4 sub-iterations which add online basis functions in the non-overlapping coarse neighborhoods u>ij with 

(i, j) £ -fa;,odd ^ -fy,odd; (b j) ^ fa;,odd ^ fy,even; (L J) £ fa:,even X Fy,odd and (i, j) £ Fa;,even ^ Iy,even respectively. 

We will take y = 2 and D = [0, l] 2 . The domain is divided into 10 x 10 uniform square coarse blocks. 
Each coarse block is then divided into 10 x 10 fine blocks consisting of uniform squares. Namely, the whole 
domain is partitioned by 100 x 100 fine grid blocks. The medium parameter k is shown in Figure [3] The 
source function / is taken as the constant 1. To compare the accuracy, we will use the following error 
quantities 

\\uh - uh\\l 2 (d) , \\uh - mf/||dg 

e 2 =-77—77-, and e a = - - — - . 

\\Uh\\L 2 (D) ||Uh||DG 
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Figure 3: Permeability field re. 

4.1 Comparison of using different number of initial basis 

In Table [T] we present the convergence history of our algorithm for using one, two, three, four initial basis 
per coarse neighborhood. Notice that, in the presentation of our results, DOF means the total number of 
basis functions used in the whole domain. We use the continuous multiscale basis functions as the initial 
partition of unity. In the tables, we obtain a fast error decay which give us a numerical solution with error 
smaller than 0.1% in two or three iterations. We can see the error decay of using one initial basis is slower 
than the error decay of using two or more initial basis since A m ; n for using one initial basis is too small. 


DOF 

€<i 

e2 

DOF 


e2 

324 

44.50% 

24.88% 

648 

17.73% 

3.58% 

648 

9.92% 

2.18% 

972 

0.31% 

1.80e-2% 

972 

0.78% 

7.54e-2% 

1296 

3.52e-3% 

1.62e-4% 

1296 

3.24e-2% 

2.13e-3% 

1620 

1.81e-5% 

8.58e-7% 

1620 

2.42e-4% 

1.10e-5% 

1948 

1.04e-7% 

4.68e-9% 

DOF 

Ca 

e2 

DOF 


e2 

972 

11.30% 

1.72% 

1296 

8.38% 

1 .00% 

1296 

0.45% 

2.44e-2% 

1620 

7.98e-2% 

3.13e-3% 

1620 

3.05e-3% 

1.37e-4% 

1944 

9.93e-4% 

3.57e-5% 

1944 

1.06e-5% 

4.08e-7% 

2268 

1.39e-5% 

5.15e-7% 

2240 

4.59e-8% 

2.14e-9% 

2540 

4.23e-8% 

1.55e-9% 


Table 1: Top-left: One initial basis (A m i n = 4.89e — 4). Top-right: Two initial basis (A m j n = 0.9504). 
Bottom-left: Three initial basis (A m i n = 1.4226). Bottom-right: Four initial basis (A m i n = 2.2045). 

To further study the importance of the initial basis, we will present another example with a different 
medium parameter re shown in Figure [4] The domain D is divided into 5x5 coarse blocks consisting of 
uniform squares. Each coarse block is then divided into 40 x 40 fine blocks also consisting of uniform squares. 
The convergence history for the use of one, two, three, four initial basis and the corresponding total number 
of degrees of freedom (DOF) are shown in Table O Table EH Table 01 Table [5] respectively. We consider 
two different contrasts. On the right table, we increase the contrast by 100 times. More precisely, the 
conductivity of inclusions and channels in Figure 2 (left figure) is multiplied by 100. In this case, the first 4 
eigenvalue that are in the regions with channels become 100 times smaller. The decrease in the eigenvalues 
will slow down the error decay. In Table O we can observe that the error decay for the lower contrast case 
is much faster than the higher contrast case. In the higher contrast case, the error stop decreasing in some 
iterations. Similar observations are obtained when we use 2 or 3 initial basis. For using four initial basis, we 
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observe a rapid convergence for both higher and lower contrast case. 



Figure 4: Permeability field k. 


DOF 

Ca 

e2 

DOF 


e2 

64 

25.44% 

6.67% 

64 

25.45% 

6.67% 

128 

1 .20% 

0.23% 

128 

1.45% 

0.27% 

192 

0.47% 

0 .10% 

192 

1.39% 

0.27% 

256 

0.26% 

5.79e-2% 

256 

0.84% 

0.15% 

320 

0 .10% 

2.30e-2% 

320 

0.34% 

7.98e-2% 

384 

6 .22e-2% 

1 .02e-2% 

384 

0.34% 

7.91e-2% 

448 

3.70e-4% 

1.57e-5% 

448 

0.15% 

3.71e-2% 


Table 2: One initial basis. Left: Lower contrast(le4)(A m i n = 0.0062). 
Right: Higher contrast(le6)(A m i n = 6.22e — 5). 


DOF 


e-2 

DOF 


e2 

128 

18.22% 

4.42% 

128 

18.56% 

4.62% 

192 

1.14% 

0 .12% 

192 

1.37% 

0.16% 

256 

0.50% 

4.95e-2% 

256 

1.25% 

0.14% 

320 

4.17e-2% 

2.06e-3% 

320 

1.23% 

0.13% 

384 

5.73e-3% 

5.38e-4% 

384 

0.41% 

3.22e-2% 

448 

7.12e-4% 

2.89e-5% 

448 

3.63e-2% 

3.56e-3% 


Table 3: Two initial basis. Left: Lower contrast(le4)(A m i n = 0.027). 
Right: Higher contrast(le6)(A m j n = 2.72e — 4). 


4.2 Setting tolerance for the residual 

In this section, we will show the performance for the online enrichment implementing it only for regions 
with a residual error bigger than a certain threshold. We consider the medium parameter shown in Figure 
El We show the results for using three different tolerances ( tol ) 10~ 3 , 10 -4 and 10 -5 . We will enrich for the 
coarse regions with residual larger than the tolerance. In Tabic El we show the errors when using 1 initial 
basis function for tolerances 10 -3 , 10~ 4 and 10 -5 . We can see that the convergence history in the first few 
iteration is similar to the result shown in previous section. Moreover, the energy error of the multiscale 
solution is in the same order of the tolerance and the error will stop decreasing even if we perform more 
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DOF 


e2 

DOF 


e2 

192 

10.69% 

1 .86% 

192 

11.55% 

2.14% 

256 

0.80% 

6 .66e-2% 

256 

1.13% 

0 .10% 

320 

0.34% 

2.24e-2% 

320 

0.98% 

8.85e-2% 

384 

1.51e-2% 

6.24e-4% 

384 

0.96% 

8.95e-2% 

448 

2.25e-4% 

1.61e-5% 

448 

0.30% 

1.39e-2% 

508 

1.72e-6% 

6.70e-8% 

508 

2.00e-3% 

8.39e-5% 


Table 4: Three initial basis. Left: Lower contrast(le4)(A m i n = 0.0371). 
Right: Higher contrast(le6)(A m j n = 3.75e — 4). 


DOF 


e2 

DOF 


e2 

248 

7.92% 

1.14% 

242 

9.63% 

1.59% 

312 

0.25% 

2.42e-2% 

306 

0.51% 

5.40e-2% 

376 

5.09e-3% 

2.72e-4% 

370 

1.38e-2% 

9.46e-4% 

440 

5.18e-5% 

2.62e-6% 

434 

2.10e-4% 

1.59e-5% 

484 

1.39e-6% 

6.40e-8% 

494 

1.74e-6% 

1.27e-7% 


Table 5: Four initial basis. Left: Lower contrast(le4)(A m i n = 0.4472). 
Right: Higher contrast(le6)(A min = 0.3844). 


iterations. Therefore, we can compute a multiscale solution with a prescribed error level by choosing a 
suitable tolerance in the adaptive algorithm. In Table [7] and Table [8l we show the errors for the last three 
iterations when using 2 and 3 initial basis functions respectively for tolerances 10~ 3 , 10 -4 and 10 -5 . We 
have the same observation that the energy errors have the same magnitude as the tolerances. 


DOF 


e2 

324 

44.50% 

24.88% 

648 

9.92% 

2.18% 

924 

0.81% 

7.72e-2% 

976 

0.29% 

2.49e-2% 


DOF 

^ a 

e2 

DOF 

&CL 

e2 

324 

44.50% 

24.88% 

324 

44.50% 

24.88% 

648 

9.92% 

2.18% 

648 

9.92% 

2.18% 

972 

0.78% 

7.54e-2% 

972 

0.78% 

7.54e-2% 

1176 

4.12e-2% 

2.88e-3% 

1284 

3.24e-2% 

2.13e-3% 

1184 

2.65e-2% 

1.57e-3% 

1364 

2.56e-3% 

1.55e-4% 


Table 6: One initial basis. Left: tol = 10 3 . Middle: tol = 10 4 . Right: tol = 10 5 . 


DOF 


e2 

DOF 


e2 

DOF 

&CL 

e2 

648 

17.73% 

3.58% 

648 

17.73% 

3.58% 

972 

0.31% 

1.80e-2% 

964 

0.33% 

1.85e-2% 

972 

0.31% 

1.80e-2% 

1248 

3.99e-3% 

1.85e-4% 

972 

0.30% 

1.63e-2% 

1136 

2.53e-2% 

1.24e-3% 

1276 

2.49e-3% 

1.19e-4% 


Table 7: Two initial basis. Left: tol = 10 3 . Middle: tol = 10 4 . Right: tol = 10 5 . 


DOF 


e2 

DOF 


e2 

DOF 


e2 

972 

11.30% 

1.72% 

972 

11.30% 

1.72% 

1296 

0.45% 

2.44e-2% 

1248 

0.50% 

2.57e-2% 

1296 

0.45% 

2.44e-2% 

1564 

3.52e-3% 

1.56e-4% 

1276 

0.24% 

9.98e-3% 

1436 

2.60e-2% 

9.70e-4% 

1576 

2.49e-3% 

1.04e-4% 


Table 8: Three initial basis. Left: tol = 10 3 . Middle: tol = 10 4 . Right: tol = 10 5 . 
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4.3 Adaptive online enrichment 

In this section, we will show the performance for the online enrichment implementing it only for regions that 
have a cumulative residual that is 9 fraction of the total residual. We consider the medium parameter shown 
in Figure 0] (4 channels medium). 

Assume that the local residuals are arranged such that 


n > r 2 > r,3 > • ■ • . 

We only add the basis </>i, • • • , (f)k for the coarse neighborhoods uq,•■ • , w*, such that k is the smallest integer 
with 

N c k 

i=1 i=l 

In Table 01 we present the error for the last 5 iterations when using 1 initial basis functions with the 
tolerance 10 -5 and 9 = 0.5. Comparing the result to the previous case, we can observe that this can use less 
number of basis functions to achieve a similar error. In Figure O we present the distribution of number of 
basis functions in coarse blocks, and see that the number of basis functions is larger near the channels (c.f. 
Figure 0]). Thus, online basis functions can be adaptively added in some regions using an error indicator. 


DOF 


e-2 

348 

0.35% 

5.19e-2% 

368 

0.27% 

4.03e-2% 

392 

6.13e-2% 

9.34e-3% 

412 

6.04e-3% 

6.60e-4% 

424 

1.51e-3% 

1.25e-4% 


Table 9: The results using cumulative errors with 9 = 0.5, tol = 10 5 and 1 initial basis. 

1 60 

50 

■ - 40 

■ - 30 

^ 20 

r 



Figure 5: Distribution of number of basis functions in coarse blocks. 


5 Conclusion 

Though the use of offline basis functions is important for multiscale finite element methods, adding online 
basis functions in some regions can improve the convergence dramatically. The construction of online basis 
functions for various applications and discretizations require a careful analysis. In particular, as we have 
shown earlier [5] for GMsFEM within continuous Galerkin framework that one needs a certain number of 
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offline basis functions in order to guarantee that the online basis functions can result to a convergence 
independent of physical parameters. In this paper, we develop an online basis procedure for GMsDGM that 
can provide a convergence independent of the contrast and small scales. Because multiscale basis functions 
are discontinuous across coarse-grid boundaries, we construct a special offline space as well as online basis 
functions. We show that our construction will guarantee a convergence independent of the contrast and small 
scales if we select a certain number of offline basis functions based on a local spectral problem. Furthermore, 
we apply an adaptive procedure to add online basis functions in only some selected regions. Numerical 
results are presented to back up our theoretical findings. 


Appendix 

In this section, we proof Lemma 12.11 Let K be a coarse grid block and let uqk be the unit outward normal 
vector on dK. We let V h (K) be the space of continuous piecewise bi-quadratic polynomials defined in K, 
and we denote V h {dK) by the restriction of the conforming space V h (K) on dK. The normal flux kVu- non 
is understood as an element in Vh(dK ) and is defined by 


/ (kVu • tiqk) ■ v = / kVu-Vw, v&V h (dK), 
J dK J K 


(30) 


/dK JK 

where v £ Vh(K) is the harmonic extension of v in K. By the Cauchy-Schwarz inequality, 


/ (kVu • hqk) ■ v < a E (u, u) 2 a E (v,v) 2 . 
J dK 


By an inverse inequality and the fact that v is the harmonic extension of v 


a%(v,v) < K K C? nv h 1 [ H 2 , 

J dK 


(31) 


where we recall that kk is the maximum of k over K and Cinv > 0 is the constant from inverse inequality. 
Thus, 


'dK 


(kVu • n dK ) ■ v < KftCinvh 2 \\v\\ L 2 ( dK) a^iu, u) 2 . 


This shows that 


/ \kVu ■ n dK \ 2 < K K C? nv h 1 o^(u,u). 

J dK 

Next, by the definition of odg, we have 

odgKu) = a H (u,v) - ^2 J ({KVtcriji}[r] + {KVo-nE}[ii]j+ ^ J k[m][4 


(32) 


Ee£t 




Notice that 


a H (u,v)+ J2 f I k H M < IMIdg IMIdg- 

EG£« hJe 


For an interior coarse edge E £ £ H , we let K + ,K £ T H be the two coarse grid blocks having the edge E. 
By the Cauchy-Schwarz inequality, we have 


I {kVu ■ n E } ■ [c] < (h J {kVm • n E } 2 {n) ^ 2 J k[vJ 2 ) 2 . 


( 33 ) 


Notice that 


h J (kVm • ?t_e} 2 (kO X <h^J (k + Vii + • n E ) 2 (n K +) 1 + J {n Vu -n E ) 2 (K K -) ^ 
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where vA = u\k±, kA = k| k±- So, summing the above over all E, we have 



< 


h V' / (kVu • HdK) 2 {nK) 1 
K£Th ^ dK 


< C? nv a H (u,u). 


Thus we have 


Similarly, we have 


y, [ {kV« • n E }[v] < C inv a H {u, u)* ( V | / k|u] 2 ds) 2 . 

\ Ee£H hJ E > 


y [ {kV-d • n B }[u] < C'i nv a_f/(u,u)2 ( y f n [u] 2 ds) 2 . 

y Ees» hjE 


Summing the above two inequalities, we have 

y / ({kVu • n B }[u] + {kVw • n B }[u]) < C inv 7 _ 5 ||u|| D g IMIdg- 

Ee£H 


This proves the continuity (f5|) . 

For the coercivity dH|), we have 


a D G ( u , u ) = ||m||dg “ E / ({«Vm ■ n E } ■ H + {kVu ■ n E } ■ [«]) • 
e<e£ hJe 


By (l35l) . we have 

which gives the desired result. 


a D c(w,w) > (1 — Cinv7 ")IMIdg. 


(34) 


(35) 
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